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Abstract 


We use the so-called t-Hill tail index estimator proposed by Fabian (2001), rather than Hill’s 
one, to derive a robust estimator for the distortion risk premium of loss. Under the second- 


order condition of regular variation, we establish its asymptotic normality. By simulation 
study, we show that this new estimator is more robust than of Necir and Meraghni (2009) 
both for small and large samples. 
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1. Introduction 


In many important applications in Finance, Actuarial Science, Hydrology, Insurance, one of 
most crucial topics is the determination of the amounts of losses of a heavy-tailed risks. In 
literature there are many possible dehnitions of risk according to the investment objectives, 
so in some sense risk itself is a subjective concept as well as the risk measure for an investor. 
Thus, the existence of a unique risk that solves the investor’s problems is not conhrmed. 
However, one must distinguish between risks. The concept of coherence due to the paper 
of (Artzner et ai, 1999) which categorize risks by good or bad risks. For this reason and 
to improve the performance of investor’s strategies we may identify those risk measures and 
the appearance of the risk himself heavy tails, asymmetries,... Most of this risk measures, 
used are special cases of Wang’s distortion premium (Wang, 1996), dehned as follows 



( 1 . 1 ) 


where '0 : [0,1] -A- [0,1] is a non-decreasing function called distortion function, such that 
"0(0) = 0 and "0(1) = 1- The distortion functions ijj are concave, which makes the correspond¬ 
ing distortion premiums H)-^; F] coherent (Artzner et al, 1999) as proved by Wirch and Hardy 
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(1999). In this paper, we suppose that the distortion functions t i —> (t is regularly 

varying at inhnity with index of regular variation p > 1, such that 

^ (t"^) = (t ), ( 1 . 2 ) 

where 1 1 —> (t) is slowly varying as inhnity, that is (tx) (t) ^ 1 as t —)■ oo, for any 
a; > 0. Note that a non negative random variable (rv) X with hnite mean and a cumulative 
distribution function (cdf) F is called heavy-tailed, if 1 — F is regularly varying with index 
— 1/7 < 0 ( notation: F G that is 

lim - - ^ for a: > 0. (1.3) 

t^oo I-F{t) ’ ^ ^ 

In particular, the proportional-hazards premium (see, Rolski et al, 1999, page 82). 

poo 

n[p;F]=/ {1-F{x)f^dx, (1.4) 

Jo 

with the concave distortion function ^p{t) = for every p > 1. Since we are concerned 
with with heavy-tailed losses with inhnite second moment, then by following Brahimi et al. 
(2011), we assume that 7 G (1/2,1) and p 7 G (0,1), thus we will work with 

1/2 < 7 < 1/p. (1.5) 

Suppose that we have an independent and identically distributed (iid) sample Xi, ...,X„ of 
rv X of size n with a cdf F satisfying condition (1.3) and let denote by Xi-n < ... < Xn:n the 
corresponding order statistics. Also, let 1 < A: = be the number of extreme observations 
used in the computation of the tail index. We assume that k satishes the conditions 

l<A:<n, A:—j-oo and k/n ^ 0 as n —)■ 00 . (1.6) 

Necir and Meraghni (2009) proposed an alternative estimator of (1.4) and establish its as¬ 
ymptotic normality by using the Weissman’s estimator of the high quantile qt = F'^ (1 ~ A) 
dehned by 

qt = {k/ny'' Xn-k:nt~'^'', 11 0, 
where denotes the generalized inverse of F and 

1 ^ ^ 

7 h (A^') • ^ ^n—i+l'.n 


(1.7) 
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is the well-known Hill estimator (Hill, 1975) of the tail index 7 . For a hxed aversion parameter 
p, their estimator is given by 

-y- n 

Hn ( 7 ^, k) := ^ ^ (1-8) 

^ P i=k+i 


The Hill estimator is a pseudo-maximum likelihood estimator based on the exponential ap¬ 
proximation of normalized log-spacings Yj = j (logXj.„ — logXj+i.„) for j = 1, In prac¬ 
tice, the Hill estimator depends on the choice of the sample fraction k and is inherently not 
very robust to large values Yj, which makes the estimator proposed by Necir and Meraghni 
(2009) sensitive. This constitutes a serious problem in terms of bias and root mean squared 
error (RMSE). To improve the quality of H^ ( 7 ^,/c), instead of Hill’s one, we propose to 
estimate the tail index 7 by the so-called t-Hill estimator, proposed by Fabian (2001), given 
by its harmonic mean 


7 = 7(/c) 





n—k:n 




-1 


- 1 


(1.9) 


known as score moment estimation (t-score or t-estimation method). The latter is more ro¬ 
bust than the classical Hill estimator 7 -^ dehned in (1.7) (see Stehlik et ai, 2012 and the as¬ 
ymptotic normality is given in Theorem 2 of Beran et ai, 2014). For other robust estimators 
for 7 we referred to Peng and Welsh (2001), Juarez and Schucany (2004), Vandewalle et ah 
(2007) and Kim and Lee (2008). The rest of the paper is organized as follows, in Section 2 
we present a construction of a robust estimator of H)-^; F] in the case of heavy-tailed losses. 
In Section 3 we establish its asymptotic normality. In Section 4 we carry out a simulation 
study to illustrate empirical performance and robustness of the estimator. Concluding notes 
are given in Section 5. Proofs are gathered in Section 6 . 

p 

Throughout the paper, we use the standard notation —)■ for the convergence in probability 
and J\f (/i, a) to denote a normal rv with mean p and variance a. 


2. Defining the estimator 

By using the generalized inverse and for a hxed distortion function '0, we may rewrite 
( 1 . 1 ) into 

n^[X] :=- [ iJ (s) dF^ (1 - s). ( 2 . 10 ) 

Jo 

The empirical estimator of the risk premium n^[X] is obtained by substituting on 
the right-hand side of equation (2.10) by its empirical counterpart F^{s) ;= inf{a; G M : 
Fn{,x) > s}, 0 < s < 1 , associated to the empirical cdf dehned on the real line, dehned 
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hy Fn{x) := n {Xj < x,l < i < n} where denote the cardinality of a set A. After 
straightforward computations, we obtain the formula 


n^X] := f 

Jo 


F^{1- s) dtlj (s) 


which may be rewritten, in terms of Xi-n, ...,Xnm, as an L-statistic 


n 


n72[^] ^ ("0) 


( 2 . 11 ) 


where 


Ci.n (i>) = <l> (tin) In). 


( 2 . 12 ) 


The form (2.11) is a linear combinations of the order statistics (see, Shorack and Wellner, 
1986, page 260). The limit behavior was discussed by many authors: Chernoff et al. (1967), 
Stigler (1974), Mason (1981), Jones and Zitikis (2003) (see its Theorem 3.2 in the case that 
X is not heavy-tailed) and in Brahimi et al. (2011) (in heavy-tailed case). 

2.1. Heavy-tailed losses case. Let X be a non-negative rv with cdf F G The 

condition (1.3) is equivalent to 



(2.13) 


= X for every a; > 0. 


We say that the function s —?■ F^ (1 — s) satisfying condition (2.13) is regularly varying at 


zero with the index (— 7 ) < 0. The parameter 7 is called the tail index or extreme value 
index (EVI). A various tail index estimators have been suggested in the literature, based for 
instance of the conventional maximum likelihood method, moment estimation, ... (see, e.g. 
Hill, 1975, Pickands, 1975, Dekkers et al., 1989. Csorgo et al., 1985 and Drees, 1995). For 
the robustness and bias reduction (see, i.e. Peng and Qi, 2004 and Stehlik et al, 2012). The 
regular-variation condition itself is not sufficient for establishing asymptotic distributions. 
To this end, we suppose that cdf F satisfy the well-known by the second-order condition of 
regular variation with second-order parameter r < 0 , that is: there exists a function t —)■ a{t) 
with constant sign at inhnity and converges to 0 as t —?■ cx) such that 



(2.14) 
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x^h — \ 

for every a; > 0. When r = 0, then the ratio-shonld be interpreted as logx. In terms 

yr 

of the qnantile fnnction condition (2.14) is eqnivalent to the following one 


lim 

ts^O 


A(t) 



T 


(2.15) 


for every x > 0, where A (t) := y^a (F^ (1 — t)), (see de Haan and Stadtmiiller, 1996 or The¬ 
orem 3.2.9 in de Haan and Ferreira, 2006, page 48). The Weissman estimator (Weissman, 
1978) of high qnantiles F^ is given by 


- s) := s i 0. 


(2.16) 


The formnla (2.10) can be split into 


nkin rl 

U^[X] = - ^jJ{s) dF^{l - s) - / ^|J{s) dF^{l - s). 

Jo Jk/n 


(2.17) 


By nsing an integration by part to the second integral yields 


= ^l:{k/n)F^{l-k/n) - 


^kln 


'i/j (s) dF^{1 — s) + / F^{1 — s)di/j (s) 

J kjn 


— -f- 


A simple estimator of H^^^ is 

•= ^n-k,n- (2.18) 

To estimate H^y we note that y is a consistent estimator for y (Stehlik et al, 2012) and since 
p < 1/y, by snbstitnting F^^’^^(l — s) given in (2.16) instead of F^(l — s) and integrating 
yield the following estimator 

pkin 

^p!n-=^i^/'^y^n-k:n (s) ds. (2.19) 

Jo 

Finally, by plngging F^ instead of F"*" on second integral of Eqnation (2.17) we obtain the 
estimator 

n 

( 2 . 20 ) 

i=k-\-l 

of where F^{s) := inf {x G M : F„ (x) > s} , 0 < s < 1, denote the sample qnantile 

fnnction associated to the empirical cdf dehned on the real line by F„ (x) := n~^ X]r=i ^ i 
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with I(-) being the indicator function and the coefficients ('^) are given in (2.12). The 
final form of our estimator 




n 




:=X„_ 


/c,n 


ij{k/n) + ^{k/ny / S ^ V (s) ds + ^ Ci,n (^/’) ^n-i+l:r 


i=k-\-l 


To establish the asymptotic normality of our estimator and compared with the estimator 
proposed by Necir and Meraghni given in Equation (1.8), we use the same function ijj (t) = 
for p > 1 used in Necir and Meraghni (2009). In this case our estimator have the 
following form 

n„ ( 7 , k) := {k/nf^'' ^ ^ (^(i/n)^/^ - ((i - 1 )( 2 . 21 ) 

i=k+i 


3. Asymptotic distribution 

We will begin to expose our results as asymptotic representations theorems in the lines of 
Beran et al. (2014). For that purpose, we need to describe the probability theory on which 
they hold. Indeed, we use the so-called Hungarian construction of Csorgo et al. (1986). 
For this we dehne {17„(s),0 < s < 1}, the uniform empirical distribution function and we 
consider the order statistics Ui-n < ... < Un-n pertaining to the independent standard uniform 
rv’s Ui, U 2 ,..., we introduce the uniform empirical quantile function {I4(s), 0 < s < 1} based 
on the n > 1 hrst observations of Ui, U 2 ,..., on (0,1) such that 

14 (s) = Ui^n for {i-l) /n < s < i/n, z = 1,..., n, and K (0) = fogn- 

Let (5n (s) = \/n (s — I 4 (s)), 0 < s < 1, be the corresponding quantile empirical process, for 
n > 1. 


We use the well-known by Gaussian approximation given in Csorgo et al. (1986) Corollary 
2.1. It says that: on the probability space (C, A, P), there exists a sequence of Brownian 
bridges {B„ (s); 0 < s < 1 } has the representation 

{Bn (s) ; 0 < S < 1} = {Wn (s) - sWn (1) ; 0 < S < 1} , 


where Wn is a standard Wiener process such that for every 0 < C < 1/2, 




sup 

l/n<s<l-l/n (s(l — s)) 


Bn(s)l ^ ( -C\ 

1/2-c - W ) 


(3.22) 


Our main result is the following 
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Theorem 3.1. Let F be a df satisfying (2.15) with 7 > 1/2 and suppose that (•) is 
continuously differentiable on [0,1). Let k = kn satisfying (1-6) such that y/kA [n/k) —)■ 0 as 
n —)■ 00 . For any 1 < p < I/ 7 , we have 

where 

2/ ^ _ 2 , 7 V (P - 2p7^ + 27 ) , _ 272 _ 

^ ^ (7P-1)^ (p + 7P-l)(p + 27P-2) 

^ 27 _ 27p (p 7^ -p7 + 1 ) 

27 - 1 (7P - 1) (P + 7P - 1)' 

4. Simulation study 

4.1. Performance and comparative study of fl^ and n„. In this simulation study we 
examine the performance of our estimator n„ ( 7 , k) given in ( 2 . 21 ) and compare it with that 
of n„ ( 7 ^, k) given in (1.8). Thus we follow the steps below. 

Step 1: We generate 1000 pseudorandom samples of size n = 100,200,500 and 1000 from 
Pareto cdf with 7 = 0.6. 


Step 2 : We estimate the tail index parameter by Hill and t-Hill estimators ^^{k\) and 
7 (^ 2 ), respectively given in (1.7) and (1.9). We adopt the Reiss and Thomas algorithm (see 
Reiss and Thomas, 2007, page 137), for choosing the optimal numbers of upper extremes ki 
and /c 2 - By this methodology, we dehne the optimal sample fraction of upper order statistics 
k* by 

1 ^ ^ 

k* := argmm-^i'^ | 7 j (i) - median { 7 ^-( 1 ),..., 7 ^-(fc)}| ,j = 1,2 

i=l 

where 71 = 7^ and 72 = 7. On the light of our simulation study, we obtained reasonable 
results by choosing 9 = 0.3. 


Step 3 : We £x the distortion parameter with respect to Condition (1.5) by p = 1.12, 
then we compute the bias and RMSE of the four estimators 7 ^(/c^), 7 (^ 3 ), and 

n„ ( 7 '^, /C 2 ). The results are summarized in Table 4.1. We see that when dealing with large 
samples our estimator performs better. 








kl 

7(/cl 

0 

( 7 , ^ 1 ) 

k*2 


^ 2 ) 

fin {l^,k*,) 

n 


bias 

RMSE 

bias 

RMSE 


bias 

RMSE 

bias 

RMSE 

100 

10 

-0.0733 

0.2511 

0.3618 

0.5199 

17 

-0.1641 

0.2865 

0.4096 

0.7332 

200 

23 

-0.0571 

0.1821 

0.3562 

0.5147 

34 

-0.0993 

0.2350 

0.3918 

0.7185 

500 

62 

-0.0299 

0.1142 

0.3404 

0.4820 

86 

-0.0301 

0.0739 

0.3639 

0.6936 

1000 

129 

-0.0147 

0.0798 

0.1966 

0.2687 

169 

-0.0181 

0.0545 

0.2827 

0.5279 


Table 4.1. n„ ( 7 , k^) and n„ ( 7 , k^) estimators based on 1000 

samples of Pareto-distribnted claim amonnts with tail index 0.6 and distortion 
parameter p = 1.12. The exact value of the premium is 2.0487. 

4.2. Comparative robustness study. In this subsection we study the sensitivity to out¬ 
liers of Yin ( 7 , k\) and compare it with that of n„ ( 7 ^, k^ . We consider an e-contaminated 
model known by mixture of Pareto distributions 

^71,72,. (a;) = 1 - (1 - e) (4.23) 

where 71,72 > 0 and 0 < e < 0.5 is the fraction of contamination. Note that for e = 0, 7 ^ 
and 7 are asymptotically unbiased. Therefore, for e > 0, the effect of contamination becomes 
immediately apparent. If 71 < 72 and e > 0, (4.23) corresponds to a Pareto distribution 
contaminated by a longer tailed distribution. For the implementation of mixtures models 
to the outliers study one refers, for instance, to (Barnett and Lewis, 1995, page 43). In this 
context, we proceed our study as follows. 

First, we consider 71 = 0.6, 72 = 2 to have the contaminated model and let p = 1.12. Then 
we consider four contamination scenarios according to e = 5%, 10%, 15%, 25%. 

For each value e, we generate 1000 samples of size n = 100, 200 and 1000 from the model 
(4.23). Finally, we compare the n„ ( 7 , k\) and n„ ( 7 -^, estimators with this true value, 
by computing for each estimator, the appropriate bias and RMSE and summarize the results 
in Table 4.2. 

As expected, the estimator n„ ( 7 ^, k^ as well as n„ ( 7 , kD turn out to be more sensitive 
to this type of contaminations. For example, in 0% contamination for n = 200, the couple 
(bias, RMSE) for If^ take the values (0.3918,0.7185), while for 15% contamination 

the bias and the RMSE are given by the couple (—1.4355,2.3107). We may conclude that 
the bias and RMSE of If^ ( 7 ^, k^) estimator are more sensitive (or note robust) to outliers. 
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Hn ( 7 , K) 

fin 

n 

% contamination 

bias 

RMSE 

bias 

RMSE 

100 

5 

0.4043 

0.6664 

-0.4286 

1.4727 


10 

0.4389 

0.6862 

-0.7291 

1.9123 


15 

0.4598 

0.7464 

-1.2786 

2.1247 


25 

1.0578 

1.1305 

-1.2103 

2.1828 

200 

5 

0.3831 

0.5532 

-0.4713 

1.5118 


10 

0.3964 

0.5675 

-1.2496 

2.1907 


15 

0.4508 

0.6870 

-1.4355 

2.3107 


25 

0.9470 

1.0197 

-1.7366 

2.3274 

1000 

5 

0.2124 

0.3211 

-0.3794 

2.1222 


10 

0.2329 

0.3349 

-1.0662 

2.3978 


15 

0.2931 

0.3749 

-1.2501 

2.0355 


25 

0.8124 

0.9291 

-1.5238 

2.3596 


Table 4.2. n„( 7 ,/c*) and n„( 7 '^,/c 2 ) are based on 1000 samples of mix¬ 
ture of Pareto distributions with tail index 0.6, e = 5%, 10%, 15%, 25% and 
distortion parameter p = 1.12. The exact value of the premium is 2.0487. 


However, for 0% contamination the (bias, RMSE) of n„ ( 7 , k^) is (0.3562, 0.5147), while for 
15% contamination is (0.4508, 0.6870). Both the bias and the RMSE of H^ ( 7 , kl) estimation 
are note sensitive to outliers. Then we may conclude that is the better estimator. 


5. Concluding notes 

We showed that the new estimator of premium based on t-Hill estimator is more robust 
and performs better than the one based on Hill estimator proposed by Necir and Meraghni 
(2009). Our estimator H^ ( 7 , k) is based on Weissman’s estimation of high quantiles, so we 
would lead to improve our result to use one of several bias-reduced estimators have been 
proposed (see for example Matthys and Beirlant, 2003). 
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6. Proofs 

To establish the asymptotic normality of 11^ „ we need the asymptotic approximation of 7 
with the same sequence of Brownian bridges as for this reason we give the following 
results. 


Proposition 6.1. Assume that the second order condition (2.15) holds with 7 > 1/2 and 
let k = kn be an integer sequence satisfying (1-6) and '/kA{n/k) —)■ 0. Then, there exists a 
sequence of Brownian bridges {B„ (s), 0 < s < 1} such that 

(7 - 7 ) = 7 (7 + 1)^ [ (s) ds + Op (1), 

Jo 

leading to 

'/kiff -^) An fo, 


Proof. Our proofs are conducted in the probability space described in Section 3. Then we 
are entitled to write 


^ 1 Xn-k-.n 




then 7 = 5'^ — 1- The asymptotic normality of 7 is established in Beran et al. 
given their Wiener process representation. Here we suppose that y/kA (n/k) A 


Ak(sk- ^ 




7 


7 + ly 7 + 1 


(1) - 7 / (t) dt + Op (1) 


(2014) by 
= 0. So 


Note that 


®n ( 5 ) = (s) - (1) , 


it follows that 

Ak (^Sk - ^ j ■ 

Using the map g (x) = Ijx — l, since g{l/ (7 + I)) = 7 and applying the delta method yields 

(7 - 7 ) = (7 + 1 )^ 7 f (s) ds + Op ( 1 ). 

Jo 

It is clear that y/k (7 — 7 ) is a Gaussian rv with mean 0 and variance . This completes 

the proof of Proposition 6.1. □ 
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6.1. Proof of Theorem 3.1. Making use of Proposition 6.1 and from Necir et al. (2007) 
we showe that under the assumptions of Theorem 3.1, there exists a sequence of Brownian 
bridges {B„ (s), 0 < s < 1 } such that, for all large n 


n(p _ 


= -7 (1 - k/n) + Op (1). 


{k/ny/pp^ (1 — k/n) 

Let U be the left-continuous inverse of 1/(1 — F). Note that U(t) is defined for f > 1. Let 

be independent and identically distributed rv’s with cdf 1 — l/y,y > 1 , and let 

bd,n < b^ 2 ,n < ••• < L^,n be the associated order statistics. Then, for p > 1, we may rewrite 
f2) 

the statistic Ll/,/^ as 


n 


( 2 ) 

ijj^n 


= {k/n)^/P-^Xn-k-,n. 

1 -7P 


Then 


1 - 7P 


where 7 is the t-Hill estimator of 7 . So we have 




where 


Ain : = 

A2n : = 
Asn : = 


{k/n)^/p\] {n/k) 

u {Yn-k-.n) 


• Ajn, 

1=1 


Vk 


IP 


1 — 7P \ U {n/k) 


y 

^ n—k:n 

n/k 


Vk- 


7P 


\/fc 


7P 

7P 


y 

^ n—k\n 

n/k 

IP 


and 


7P 1 - 7P 

(k/n)^/P p 


A4n := \fk{k/n) 


g7^U(n/fc)-n5( 


( 2 ) 


U {n/k) 

As showed in Necir et al. (2007), we have : Ain —t 0 and A 4 n —?■ 0 as n —)> 00 . 

Next, we show that A 2 n +A 3 n is asymptotically normal. Assume, without loss of generality, 
that the rv’s (yn)n>iare defined on a probability space (f 2 . A, P) which carries the sequence 


{Un)n>i in such a way that = {1 - Un) ^ for n = 1,2,... and Fj^n = (1 - ^i,n) \ i = 

l,...,n,. Then, this allows us to write Fn-i+in = (1 — fn (1 — s))~^, for - < s < —, 

’ n n 

i = 1, ...,n. From Necir et al. (2007) we have For A 2 n 

2 

A2n = - {n/ky^"^ — (1 Op (l))Bn (1 - k/n). 

1 -7P 


-1 
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For A 3 „ and by using the map h{6) = pj — p) and applying the delta method yields: 


^3n — 


7- 

(P7 - 1) 


From Corollary 6.1 we get 


'^3n — 


IP (7 + 1)^ 

(P7 - 1)^ 



s'^' (s) ds + Op ( 1 ). 


Finally we have 


(ng!. - ng') 

{kjnYIPF*^ (1 — /c/n) 


^ (s) (is-(n//c)^^^ ——B„ (1 -/c/n)+Op (1). 

(P7 - 1) Jo 1 - 7P 


From Necir et al. (2007) we have 

(<!. - ng*) jy, (l-s)F^'(l-s)ds 

(klnylPF^ (1 - kin) pp^ (i _ k/n) {klnp'‘~''^ ^ 


Then 


\/^ fn^ Yi-ijjj 

{k/ny/p-F^F^il -k/n)~^ ^ 

where 


^ (7) p) •= WAl + hFn2 + hF„3 + Op (1) 

and 


:= - {n/kY^^'yMnil - k/n) , 


W^2 



CO 

pF^ (1 — k/n) {k/nY^^~^^'^ 


It is clear that A ( 7 , p) is a Gaussian rv with mean 0 and variance 


B(A(7,p)f = E(wP)+E(Wl,^+E(Wl,^+2E(W„iWni) 

+2E (WmW„3) + 2E (Wn2W„3). 
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An elementary calculation gives, we get 


E (W^i) 

E {W^2) 


7^ + o(l), 




E{Wr,lWn3) 


E{WnlWn2) 



and 




(7P- 1) (p + 7p- 1) 


The proof of Theorem 3.1 is completed by combining all the preceding results. 


□ 
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